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Abstract: Based on quantum origin of the universe, in this article we find that the universal wave 
function can be far richer than the superposition of many classical worlds studied by Everett. By 
analyzing the more general universal wave function and its unitary evolutions, we find that on small 
scale we can obtain Newton’s law of universal gravity, while on the scale of galaxies we naturally 
derive gravitational effects corresponding to dark matter, without modifying any physical principles 
or hypothesizing the existence of new elementary particles. We find that an auxiliary function having 
formal symmetry is very useful to predict the evolution of the classical information in the universal 
wave function. 
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1. Introduction 


The principal methodology of quantum physics is to find appropriate wave functions 
for complicated quantum systems and manifest their fundamental properties. Finding ap- 
propriate single-particle wave functions can greatly enrich our understanding of complicated 
systems, such as the eigenfunctions of infinite potential well and periodic potential for studies 
of crystal. Moreover, finding appropriate multi-body wave functions proved to be helpful for 
our understanding and design of complicated quantum material; for examples, the BCS wave 
function in superconductor [1] and Laughlin wave function [2] in fractional quantum hall 
effect. The most complicated wave function we can imagine ought to be the universal wave 
function. Everett analyzed formal solutions for the universal wave function [3] and presented 
predictions which changed our understanding of the universe, such as the predictions that 
there are many classical worlds in the universe and those classical worlds are constantly 
branching. Unfortunately, the formal solution given by Everett does not have observational 
and falsifiable predictions and thus it failed to become a realistic physical theory at the present 
stage. 

In this article, we think that the universal wave function can be much richer than Everett's 
formal solution, by taking adequate consideration of quantum origin [4] of the universe. We 
discover that the diffusion particle wave packets due to quantum origin of the universe can 
have gravitational interactions with our classical world. We propose a general method to 
extract quantum information in the universal wave function; by using it we can recover 
Newton’s law of universal gravity on small scale, while on the scale of galaxies we obtain the 
gravitational effects corresponding to dark matter, without modifying any physical principles 
or hypothesizing the existence of new elementary particles. We also discover the usual 
Newtonian gravitational constant can be different from the gravitational constant between 
elementary particles. Our work clearly shows that the universal wave function is not only 
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realistic, but we can also obtain wave packets distribution information of the universal wave 
function through gravitational effects. 


2. The Unitary Evolutions of Universal Wave Function 


Mathematically, if we apply quantum mechanics to the entire universe, then at least 
under non-relativistic framework, we can assume the existence of a universal wave function 
Yu(f1,--:,?n,t) [3]. Since we are working under non-relativistic framework, we do not 
consider those zero-mass particles, nor do we consider the creation and annihilation of 
particles. The universal wave function satisfies the linear unitary evolution of quantum 


mechanics: 
Ov» 


ih 5 


where Ho denotes the kinetic energy operator of all the particles, Avsw denotes interactions 
except gravitational interactions. H. ¢ denotes gravitational interactions. 





(71,- 7 TN, t) = (Ho =p esis + Ag) ¥u(Fi,- ia TN, a (1) 


ai =, N mmMn 
g=-5 


feta Fi Fal . 
where G denotes the gravitational constant between elementary particles, m) is mass of the [th 
particle. Another equivalent way to consider the above Schrédinger equation is the Feynman 
path integrals for the universal wave function [5]. 

Even if we ignore the internal freedom of elementary particles, the universal wave 
function will still be extremely complicated. Therefore we need to use suitable mathematical 
tools to analyze the universal wave function. Mathematically, we can define the single-particle 
wave packet distribution function as 


it) =f TL Put FOR 3) 
i(Aj) 


p;(7;,t) reflects the wave packet distribution of the jth particle in the universe. Larger p;(7;, t) 
implies that more wave packets are distributed at that location. Naturally we can generalize 
the above definition to two-particle joint wave packet distribution function as 


Pin (T1Th,t) = I] @7il¥u@--- Fy, t)/?- (4) 
i(ALh) 


One may also generalize the definition to general n-particle joint wave packet distribution 
function (}.... n(71,°°* ,Tn,t). 


3. Classical World Wave Function, Classical Information and the Universal 
Wave Function 

To analyze precisely universal wave functions containing the classical world, we must 
first understand the fundamental properties of wave functions behind the classical world. In 
our classical world most particles are formed by quantum correlation with nearby particles. 
Here quantum correlation means that the wave function of two or more particles cannot be 
written in a product form. For example, in hydrogen atom, (Fp, 7e,t) cannot be expressed as 
(7p, Fert) = Pp (Fp, t)he(Fe, t), this implies there is quantum correlation between the electron 
and proton within a hydrogen atom; the same holds for the three quarks within a proton, 
which is why we regard hydrogen atoms and protons as composite particles. Given that the 
quantum correlation between the three quarks of a proton is extremely robust, we know the 
rest mass of every proton is the same and can be regarded as a fundamental physical constant. 
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There should be a wave function to describe a given classical world in which particles 
are distributed in localized regions. We call such a wave function the classical world wave 
function ¥.(7,--- ,7N,t), to distinguish it from the universal wave function ¥,(71,--- ,7N,t). 
From our understanding of the classical world, the classical world wave function should 
satisfy the following necessary conditions. 


1. The first necessary condition satisfied by ¥¢(r1,--- ,/n,t) is: almost every particle in 
Y¥c(71,--: ,7n,t) has some sort of quantum correlation with one or more other particles. 
Such quantum correlations distribute over the entire universe in a network form in our 
classical world. 

2. The second necessary condition for ¥¢(71,--- ,7Nn,t) goes like this: in our classical world, 

almost every particle is localized and distributed in a small region ©;(t). 
For classical world, we require that the volume of &;(t) be on microscopic scale for 
most particles. Under this condition, the sets {2;(f),j € {1,---,N}} and {7((t),j € 
{1,---,N}} represent the classical information in the classical world within the universal 
wave function, here 7;(t) is about the center of region &;(t) and defines the classical 
location of a particle. The gravitational effects we observe in reality correspond to 
evolution of the set {7;(t),j € {1,--- ,N}} over time. 

3. Finally, the third condition for ¥¢(71,--- ,7N,t) is as follows: when 7; and 7, are separated 
by macroscopic distance, the quantum correlation between those two particles can be 
neglected. More precisely, we have the decomposition 0) , ~ 01n- 


From the above three necessary conditions, we can see that the classical world wave 
function can only reflect partial information of the universal wave function, since in the lo- 
calized particle regions ¥.(71,--- ,7n,t) = Yul? € X(t), --- Fi E Zy(t), ++, Fn € Un(f),f), 
whereas in other regions ¥-(71,--- ,7,t) = 0 based on our request of the classical world. This 
shows clearly that ¥.(71,--- ,?n,t) A VYulf1,:-- ,7n_t). Itis clear that ¥,, is more fundamental 
than ¥.. 

Generally speaking, the universe wave function is much more complex than the classical 
world wave function, considering the quantum origin [4] of our universe. Hence, we need 
a general way to get the classical world wave function and the corresponding classical 
information if we know the exact universal wave function. Here we present a general method 
to search for classical worlds of a given universal wave function. 

Choose two arbitrary particles i and j, compute 9;(7j,t) and ;(7;,t) from 
Wul(?1,-:+,?n,t), where we assume 9;(7;,t) and Pj (7, t) each has P(t) peaks. Now we can 
construct P?(t) different functions as follows: 


cui(t)¥ul?is-++ Fut) — if F © UE(t) and 7; € Xi(t) 


spike chee — (6) 
0 if 7, € Uj (t) or Fj € Xj(t) 


Fourey (Fi, 0+ + 7 t) = 


Here cx,(f) are normalization constants. From the above equation we have constructed a 

set of functions with P?(t) elements: tigen (Tit, na ye  & qe PO), 
mj 

te {1,---,P(t)}}. 

For a given (x,1), if ¥pxyy)(71,--- ,?n,t) has about N peaks in {2,(t),j € {1,---,N}}, 

mj 

and satisfies the three necessary conditions for the existence of classical worlds when all 7; 
are confined to the domains %j(t), then it is possible that we have found a classical world. 
Additional conditions may be imposed to check if that is indeed a stable classical world 
wave function. After we get a classical world wave function, the classical information 
{21 (t),--- ,Un(t)} will be obtained simultaneously. 


Symmetry 2021, 13, 193 


4 of 18 





Of course, the above extraction of the classical world wave function and the correspond- 
ing classical information based on two particles can be generalized to n particles if two 
particles are not enough to get the definite classical world wave function. 

After the definition of the classical information based on the universal wave function, the 
question is how we determine the evolutions of the classical information. At first sight, we 
would adopt the following “type-I” program. 


1. At fo, we get the classical information of the classical world from the universal wave 
function ¥,, (to) with the above method to extract the classical information. 

2. To predict the evolution of the classical information, we use Newtonian mechanics to 
calculate the evolution of the classical information {7;(t),j7 € {1,---,N}}- 


Since there is no external observer to interact with particles in the universe and record 
their classical information, we should always analyze the universal wave function from 
unitary evolutions. When the universal wave function is applied to the whole universe, we 
should not use the concept of instantaneous wave packet collapse. It is clear that to consider 
in the most accurate way the evolution of the classical information, we should first solve the 
unitary evolution of the universal wave function, and then search the classical world from the 
universal wave function at different times. Hence, the improved “type-II” program should be: 


1. At fo, we get the classical information {7;(to),j € {1,---,N}} of the universal wave 
function ¥, (to) by searching the classical world. 

2. From the solution of the Schrédinger equation, we get the universal wave function ¥ , (tf) 
at a later time. 

3. From ¥,,(t), we get the classical information {7;(t),j € {1,--- ,N}} by another search 
of the classical world. 

4. Summarize the classical law from the relation between {7;(to),j € {1,---,N}} and 
{Fi(t),j € {y+ NP}. 

Of course, the “type-II” program will provide more accurate evolution of the classical 
information, compared with the “type-I” program. However, at first sight, it seems that there 
should not be significant difference between the “type-I” program and “type-II” programs. 
Nevertheless, we will show that our universe provides a natural example to demonstrate the 
difference between the “type-I” and “type-II” programs. 

Since we do not know what the universal wave function looks like and how it evolves, 
the “type-II” program is not as useful in practice without special mathematical technique. In 
the following, we will attempt to construct a function ¥(f), which we call the x function, 
which unlike ’,(t), takes into account the evolutions of the universal wave function itself and 
provides better approximation to predict the evolution of the classical information we observe 
in reality, compared with evolutions of the classical world wave function by Schrédinger 
equation alone. This ¥,(f) function has the merit that it can naturally take into account the 
observational data to predict the details of the future evolution. 


4, An Analytical Approach to the Universal Wave Function 


In the 1950s, Everett [3] analyzed the following formal solution for the universal wave 


function: 
P(t) 


Wulfi-+: ,Tn,t) = So a(t) PE (Fi, ,?n_t). (6) 

K=1 
Here every (tf) represents the classical world wave function and Equation (6) expresses 
that there are P(t) different classical worlds. Normally, people think the superposition of 
classical worlds constitutes the whole universal wave function. However, we consider that 
to be a big mistake. The viable formal solutions of the universal wave function can be much 
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richer than Equation (6). We will introduce a general method to give more general universal 
wave function in due course. 

For the formal solutions given by Equation (6), it is sufficient to consider the evolutions 
of ‘¥-(t) alone in Schrédinger equation, and we will prove in Section 5 that the “type-I” and 
“type-II” programs will give the same evolution of the classical information. However, in the 
following we give a simple and realistic example to show that in general this is not the case. For 
simplicity, we consider the case of a single classical world wave function. Figure 1 illustrates 
p; of the universal wave function for several particles. There is a common misunderstanding 
that ¥,, = 0 for any 7; ¢ Xj(t). Everett’s formal solution (6) clearly satisfies this property. 
However, considering quantum origin of the universe, we cannot exclude the possibility that 
p; also has small background distribution. In this case, it is possible that '¥,, 4 0 for 7; € 4;(f) 
and all other 7, ¢ X;(t). So in this case, we should consider to extend the domain of Y,, to the 
entire universe to consider the long-range gravitational effect in a more accurate way. 


P2 P3 P4 Ps 





2 3 4 5 


Figure 1. As an example to illustrate the universal wave function and its classical world, we show the single-particle wave 


packet distribution for five particles. The quantum correlation between adjacent particles is also shown. We emphasize the 


possibility that there is a nonzero background distribution outside the region of the classical world; we can imagine this by 
thinking that at first the universal wave function is the product of the diffusion wave packets of those particles. Then we turn 
on interactions to form the classical world. It is clear that without an external observer, there must be a significant amount of 


diffusion wave packets in the universal wave function. The main result presented in this work is that those diffusion wave 


packets will influence the classical world through gravitational forces. 


As an example, we first consider the domain extension of the jth particle and define the 
function ¥}(71,--+ ,Fy,t) so that 7) € )(t) for all! # j, whereas there is no confinement on 
7;. In this case, we may rewrite the normalized That) =e (t)$;(7;,t). 
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Here ey (t) is defined as follows: if Fj is inside localized region 4; (t), then ey! (71,+-+ , Fx, t) 
= ¥.(r1,--- ,7n,t), whereas if it is outside region L(t), then PLS . ile ++ Py, t) = 


Lo(P1, Fi (t), Pn, t)/ y/p; (F(t), t). For 7; € X;(t), note that 7;(t) denotes center of region 


Xj (t) rather than the independent variable 7;. For 7; ¢ £;(t), the number of independent 


variables in ¥).(7,--- /Fi-1/7i,+++ Tn, t) is3(N —1) +1. p;(7;(£), £) is the single-particle wave 
packet distribution function of '¥.(t) at location 7;(t). 

$;(7,t) has the following properties: (i) when |7; — 7;(t)| < Rj, |Gj(#,t)/? = y(t); Gi) 
when 7; — F;(t)| > Rj, bi (7;, th= (Gilt) (7, t), where pj (Fj, t) is normalized and diffuses 
across the entire universe. In the latter case, |@;(7;, t)|? is very small, so we can freely choose 
Rj on macroscopic scale, for example one meter, as long as it is greater than the radius of the 
localized wave packet distribution region 4;(t). Moreover, since |$;(7;, t) | is very small and 
R; is of macroscopic scale, |¢;(7;, t)|* is not very sensitive to the exact location of 7;(t). The 
normalization condition requires that y;(t) + ¢j(t) = 1. 

We may consider the domain extension of other particles in the same way. It is clear that 
there is no quantum correlation between ¢;(7;,t) for different particles. Proceeding in this 
way, we obtain the following form for Y(t) by the domain extension of all 7; based on ¥¢(t): 


Ex(Fi- ++ Ft) = ¥e(t)hr (71, t) G2 (72, t)...P; (Fj, £)...PN (Fn t), (7) 


where all the $j (Fj, t) are defined in the same way as above. For ¥.(t) in this case, however, 
we need to consider for each particle whether it is in its localized region %j(t). That is, we 
need to define ¥,(t) as follows: if all the particles are inside their localized regions &;(f), 
then ¥,(7,,-+-,7y,t) = ¥c(71,-++ ,Fn,t); ifa subset A (say A = {1,2,...,n}) of particles are 
outside their localized regions & j (t), while all others are inside their localized regions, then 
Ve Fnt1, isp? Ny t) =" (71 (ee Tn (t), Fats sig EN t) / V Pi,..n (1 (2), Tn (t), ay where P1,...n (t) 
is the joint wave packet distribution function of ¥,(t) for those n particles. We can obtain 
similar formulas for each subset of particles and together they constitute complete definition 
of ¥,(t). 

It may seem that the x function ¥,(t) was transformed from ¥,(t) and ¥,,(f) by a 
rather artificial mathematical operation, in the following we elucidate why the x function 
we introduce here has physical origin. In the early universe the wave packet of each particle 
completely spread out across the universe and there was no definite stable structure formed, 
no classical world existed at that point. Then as the universe expands, wave packets continue 
to spread and from the unitary evolutions of universe, we know those wave packets cannot 
“collapse” based on unitary evolution, which means that the diffusion wave packets always 
exist in the universe. So how did the classical world take its form? Because of the short and 
long range interactions between those particles, quantum correlation will stimulate structures 
to be formed among the diffusion wave packets and the network structure of quantum 
correlation between nearby particles represents the classical world. In general p;(7j, t) for 
each particle looks like a series of spikes with a diffusion background, each spike represents 
a distinct classical world. So as long as the diffusion background is not equal to 0, we can 
imagine the existence of Y(t) by extending the domain of ¥-(t) to the entire universe. 

From Equation (7), we see that information about ¥;(t) is contained within Y(t) with 
no changes. However, ¥,(t) contains more information through {7;(t),¢j(t), Pj(7j,t),7 € 
{1,---,N}}. In this work we consider a random model in which all ¢;(t) are random numbers 
(noise) between 0 and 1, while {;(7;,t),7 € {1,---,N}} are a series of functions with no 
correlation with each other as time proceeds. Considering the complexity of the universal 
wave function and its quantum origin, and considering the rapid movement of any particle in 
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the classical world relative to the cosmos background, at each time we need to reconstruct ¥y 
from ¥, and ¥., which makes this random model reasonable. Of course, in this model, there is 
no unitary evolution relation between Y, at different times. Hence, there are only continuous 
evolutions of a single classical world (if we do not consider the classical world branching 
of ¥, by itself), even when ¥, # ¥-. The correlation time of y;(t) can be estimated by 
considering that when there is no overlapping between &)(f) and X;(t+ T) (/ denotes a nearby 
particle), {y;(t), @j(7,t)} and {7j(t + T), oj(7j,t + T)} may be assumed to be independent. 
Based on this assumption, it is estimated that tT ~ 10~}° for a particle in the Milky Way. In 
this estimation, we assume that the velocity of the Milky Way relative to cosmos background 
is about 200 km/s, while the size of £)(t) is about 10~!° m. 

If at time t, ¥,(t) = Y(t), then it is not the random model we consider. However, 
we believe that the universal wave function is much richer than this simple case so that ¥, 
may be significantly different from Y, and ¥,, at any time. Of course, only observation will 
determine which model is correct, because it is impossible to solve the Schrédinger equation 
for the universal wave function to prove this. We will show in due course the observational 
evidence of our model. 

Notice that, in our model ¥y is not a real wave function at all; it only provides a mathe- 
matical method to consider the correlation between our classical world and other parts of the 
universal wave function, which enables us to study its gravitational effects later. It is nonsense 
to calculate the long-time unitary evolutions of ¥,(f). Only the unitary evolution of ¥,,(f) 
is always correct, while the unitary evolution of ¥.(f) may provide approximation to true 
evolutions of the classical information. We will show that Y(t) gives better predictions to the 
gravitational acceleration of the particles at time f in the classical world than those obtained 
from the unitary evolution of ¥,(t). 

Here, we show two more properties of ¥ ,(f). From the general expression of Y(t) given 
by Equation (7), we have 

(1) For normalized ¥,(t), 0;(7,t) = Y(t) (Fj, t)+ Cj (t)p4 (7;,t) for the jth particle. Here 
the normalized pj (7j,t) denotes the localized part distributing over region ¥j(t), whereas 
p% (7;,t) is the diffusion part. 

(2) For the original ¥.(t), when two particles are separated by macroscopic distance, 
Pim(T Fmt) = pF (71) On (7m,t). For ¥ y(t), there can be no quantum correlation between 
diffusion wave packets of different particles, so we can assume that the decomposition formula 
still holds, namely 


Pin(Fi Fn, t) = (qa(t)pi (Ft) + Sr(E) pF Ft) (Yn(E Punt) + Sn(t)en(Fnt))- (8) 


Before considering gravitational interaction, we first think about whether the existence 
of ¥, (A ¥c) will influence evolutions of the classical information in ¥;(f) through Hesw. In 
the classical approximation, evolutions of the classical information of a particle is manifested 
by acceleration 7 (t). The most useful and efficient way is to calculate < | Bess |¥, > first. 
Since the diffusion wave packets diffuse very sparingly over the universe, they will not cause 
short range interactions such as weak and strong interactions, and because the diffusion wave 
packets as a whole is electrically neutral, they will not cause electromagnetic interactions either. 
That is, < ¥y(t)|Hesw|¥ y(t) > will not contain the information 1Oltat) fe tle Nf in 
¥,(t). This leads to the result that the diffusion wave packets will not influence evolutions of 
the classical information in Y(t) through Hesw. We will show in the following section that 
gravitational force is different. 
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5. Gravitational Acceleration of Particle Wave Packets in the Classical World 


To obtain the acceleration of a particle by inter-particle interactions in the classical 
approximation, the most useful and efficient way is to calculate the interaction energy first. 
Quantum mechanics tells us that, for a given state Y(t), the total gravitational interaction 
energy is determined by the joint wave packet distribution functions as follows: 


y pet 
E,(t) = (¥(t)|Ag|¥ (8) )=-56 Damm, [and n ae (9) 
lén n 


From the above formula, it is straightforward to obtain gravitational effects in the classical 
approximation. Here 0) ,(71,7n,t) is the two-particle joint wave packet distribution function 
of ¥(t), while A, is given by Equation (2). 

If we calculate E,(t) of ¥,(t), it is clear that E,(t) will contain all the information 
{7 (t), 6) (£), Gi (7), t),7 © {1,--- /N}} in ¥,(t), because the gravitational interaction is univer- 
sal belwieen any two particles. Hence, it is expected that gravitational acceleration 7;(t) in the 
classical world will be influenced by the diffusion wave packets contained in ¥ , (tf). We will 
use the usual method to calculate the gravitational acceleration. For the wave packet of the 
jth particle in 4; (+), we consider an infinitesimal displacement A7; of this wave packet with a 
group accolanation Z a;(t) during an infinitesimal time interval Ab Based on the fact that the 
increase of the average kinetic energy is the decrease of the gravitational interaction energy, 
we can get the expression of the group acceleration dj(f). 

First we consider the simplest case: suppose the universal wave function consists of a 
single classical world and there are no diffusion parts, ie., ¥,,(f) = ¥-(t) = ¥,(t). In this 
situation, when two particles are separated by macroscopic distance, we have 


Pin(F1, Pn t) = pr(F1,t) Pn (Fn, t). (10) 


Substituting into Equation (9), we get 





= mim 
Ps) ~ ~9 OR HOT = 


Now we analyze the gravitational force for the jth particle. We assume that during an 
infinitesimal time interval At all of the momentum changes for this particle are caused by 
gravity. Assuming that the wave packet of this particle moved from 7;(t) to 7;(t + At), while 
its group velocity changed from d;(t) to dj(t + At). Then, based on energy sonseryetion, we 
have 


= jee ae - = 
m; (|Bj(t + At)|" — |8j(t)?) = Bg (F(t), ++) — Eg F(t + At), -). (12) 
First we consider the x-components oF and ir the left hand side of the above equation 
becomes: 
1 x 2 40; 
5mj( (f(t + At)? — (of (t))?) = mjot tat. (13) 


Whereas the right hand side is: 


7 7 dE 
Eg(F;(t),-++) — Eg (F(t + At),-+-) =— are (14) 
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where Ar; = r (t + At) — : (tf). Since oF - Ar; / At, we have the group acceleration 


x 
_ 40; 1 «Es 
] dt mj or; . 








(15) 


Using Equation (11), we obtain the following acceleration formula in the classical approx- 
imation: 





> mi; = = 
1 be Sy — FRO HO (16) 
This is the equation of Newton’s law of universal gravity. 
Now we analyze the formal solution (6) given by Everett. In this formal solution, every 
Yt represents a classical world and there are still no diffusion terms, i.e., if we perform the 
domain extension for any ¥¢(t), we have ¥¢(t) = ¥}(t). The two-particle joint wave packet 
distribution function of ¥,,(£) becomes: 


P(t) 
Pigltitel) => alt) ora or Aut), (17) 


K=1 


where normalized p/(7),t) represents the single-particle wave packet distribution function 
for the /th particle in the xth classical world within region LF (t). Thus, the total gravitational 


interaction energy of Y(t) is 


eo(t)~ 25 la(oe(Eoe et _). 18 
g(t) 2 Dy lal )| ve [F*(t) — 7*(t)| ( ) 


ifj i j 


Note that there are no gravitational interactions between different classical worlds in this 
case, because there is no gravitational interaction energy between different classical worlds. 
For this universal wave function, we can obtain the gravitational acceleration formula 
for the jth particle in the «th classical world from the following relation: 
1 2 aK 2 
5 la (8) mi ( a (t+ At)| = 








(0) = Ey(F(t),-++) — Eg (E+ AN,---). (19) 


|x,.(t)|* on the left and right hand sides cancel out, so we have 





t= Yi Cam — appt) — FO). ee 


This implies each classical world evolves independently according to Newton’s law 
of universal gravity and such evolutions do not depend on |a,(t)|*. Based on Everett’s 
universal wave function (6), at macroscopic scale, we can not notice different gravitational 
effect, compared with Newton’s universal law of gravity. 

For the following uncorrelated diffusion universal wave function, 


P(t) N 
Pult) = Dae (t) Pe (71,-++ Fn, t) + Blt) [1 ¢7,), (21) 
K=1 j=l 


similar calculations show that the gravitational accelerations are still given by Equation (20), 
since the gravitational interaction between each classical world and f(t) Te 1 CH (7, t) is 0. In 


the above expression, we assume that all CH are diffused in the entire universe. 
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In Everett’s universal wave function given by Equation (6), for each ¥*(71,--- ,7,t), the 
domain is highly localized. We have emphasized previously that this is a very specialized 
universal wave function. When diffusion wave packets are fully taken into account, we 
should use ¥',(f) to calculate the gravitational acceleration at t. For a given ¥',(t), based on 
the formula given in this article, for two particles separated by macroscopic distance, we have 
the two particle joint wave packet distribution function given by Equation (8). 

For the case of a single classical world, we may always decompose ¥,,(f) as 


Yu(t) = ay(t)¥ x(t) + BE) Fr (Et). (22) 


Here the term f(t)¥;(£) is due to the general consideration that ¥,,(t) A ¥,(t). For any 
Y, and Y,, we can always make this decomposition. We may assume further that 


< ¥,(t)|¥,(t) >=0. (23) 


It is clear that the superposition principle of quantum mechanics can always make the 
appropriate decomposition (22) to satisfy this request. 
The total gravitational interaction energy of ¥,,(t) is therefore: 


Eg(t) * lay (#)[?E (t) + |BCPES(t). (24) 


Here ee (t) is the total gravitational interaction energy of ¥,(t), while E¢(t) is the total 
gravitational interaction energy of ¥,(t). It is clear that the term < ¥y(t)|Hg|¥;(t) > is 
negligible. The above expression of E,(t) is the motivation of the introduction of the x function 
Y to calculate the gravitational effect in the universal wave function. From Equations (8) and 
(9), we get 





5 howl mm; 
aes) Yi iat —FID 
7 Ome) 
su jt) 
~Lewi(tGjlt ’ err 05 
ifj i Iri(t Fi — 
> d/z 
Be, mipf (Fi, t) mio! (7), £) 
- 3 bon) g(t) fa? 7 | IFi — 7; ‘ 


where 7;(t) denotes the classical position of the jth particle within region 4; (t), whereas 7; 
deriotes: position variable to be integrated. 

Now we consider the gravitational force for the wave packet of the jth particle in the 
classical world, i.e., in the region & j (t). Again, we assume during infinitesimal time interval 
At, all momentum changes for this wave packet are caused by gravity, we have 


alt) 5m (loi(t +0)? — |8)(8)7) = ERG (#),-+-) — ERG EF AL, +). 26) 


For oF and Vr we have 





do 
ay(t)gm;((of(-+ A8))? — (@F(O)?) = 4(t)mjoe hat 27) 
and 
X(F x OES x 
EX (F(t),---) — ERG; (t+ A8),---) = — Sart. (28) 
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In this case, we have 














do; 1 dex 
} J] _ g 
8 ae > ag fm ort (29) 
Therefore, 
= i(t)m; 7 ill mio Fit) 2 
= = oO a yp AO ao 6 fea FHP (Fi —F(t)). (30) 


It is worth pointing out that the value of 7;() does not appear in the gravitational 
acceleration d@j(t) of the jth particle. Hence, the fluctuation of the gravitational acceleration of 
a particle is due to the fluctuations of ;(t) of other particles. 

We assume that the probability distribution of the value of ;(t) at different times is the 
same for every particle. Hence, we may use the time average values 7 and ¢ to replace 7;() 
and ¢;(t) in the above equation. 7 and ¢ can also be regarded as the averages for relevant 
particles. When gravitational effects due to the fluctuations of 7;(t) and ¢;(t) are omitted, we 
have 


_, dmipi(? (Fa t) > > 
+ Gn f Rap) 
i(#j) 





i= Le Gua hl -F 


i(Fi) (31) 


where the second term denotes the gravitational acceleration due to diffusion wave packets in 
Y(t). In the above equation, r 


(d@+1)° 

The parameter d = ¢/7 characterizes the gravitational interaction due to the diffusion 
wave packets in ¥y(f). 

For particles in the classical world, from a gravitational viewpoint, the second term in 
Equation (31) suggests that the equivalent mass of the diffusion part is thus My = dM,,, where 
My, = a mj is the sum of the actual masses of all the N particles in the universe. Since the only 
way in which the diffusion wave packets can affect the classical world is through gravitational 
interactions, we call M, "fictitious" gravitational matter. 

When we observe gravitational effects in a small region, the contribution from diffusion 
part can be neglected and the gravitational constant we obtain through such observation is 
Gy, which we call Newtonian gravitational constant. On the other hand, the constant G is a 
constant related to elementary particles and is different from Gy, we call G the gravitational 
constant of elementary particles. After we obtain Gy, we can derive G through G = (d+1)Gy. 


Gn =Gy= (32) 


6. Applications to the Gravitational Effect of Dark Matter 


The mathematical model of ¥, provides a way to determine the size of d. Therefore, d 
becomes a parameter originating from many-body evolutions of quantum mechanics, like 
lattice constant of a crystal. We cannot derive d from first principles, but we can experimentally 
measure the size of d. 

The method for measuring d is as follows: first based on Equation (31), we can omit 
the second term and measure the value of Gy from earth-based small macroscopic scale 
experiments. Next we analyze the motion of objects on a larger scale, such as the motion of 
objects on the edge of galaxy. Therefore we can obtain d and Mj, from the difference between 
the observed a; and the theoretical result of the first term in Equation (31). Fortunately, the 


J 
remarkable advances on the astronomical observations make this scheme completely feasible. 
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From the value of Gy and d by observations, we get the true gravitational constant G for 
elementary particles. 

As we all know, gravitational effects for objects on the edge of galaxy are not determined 
by the first term in Equation (31) alone [6]. In fact, a; calculated from the first term in Equation 
(31) is several times smaller than the observed d;. Up to now, people still could not explain 
this difference and they hypothesized so called dark matter. Now we naturally discover the 
possibility that, the “gravitational matter” discussed here can be fit into the dark matter model. 

If we define the classical world mass density distribution function by 


pm (Ft) = Yomj% (Ft) (33) 
J 


and gravitational matter mass density distribution function by 


pa(7,t) = d) | mpi (7, t), (34) 
i 


then based on all the observational and theoretical researches of dark matter [6], we can 
reasonably hypothesize that o,4(7,f) may have large scale structure. In addition, it seems 
that we may regard p,(7,t)/d as the mass density distribution of Y(t). This means that 
from the observed “dark matter” distribution o4(7,t) and the numerical simulation of the 
evolution of the universal wave function, we have the chance to know the initial wave function 
of the universe, which plays a special role in the initial state [7] of inflation theory [8]. Of 
course, to give an accurate deduction of the initial condition at the inflation stage, we need 
to improve our theory to the relativistic version of quantum gravity, such as the possible 
application of the Wheeler-DeWitt equation [9]. It is worth pointing out that the universe 
wave function considered in this paper is different from that of the wave function of the 
universe in [10], where the formal solution of the Wheeler-DeWitt equation is considered 
to reveal quantum gravity. In our work, however, we do not consider the quantization of 
gravity; what we consider is the gravitational effect due to the diffusion wave packets after 
the quantum origin of our universe. Hence, in Equation (31), the right hand side does not 
contain a term comprising fi term, while we expect that further improvement by including the 
quantization of gravity will lead to a correction of f term. 

From the present observational results, we have d ® 5.5 [6], this means the diffusion 
part can be rather significant. Under non-relativistic framework, the total particle number 
N of the universe is conserved. Thus in a given classical world, no matter how we count 
the number of total particles with a detector, the detectable total particle number is always 
N. Therefore, for observers in the classical world, the total mass of the universe they obtain 
by counting the particle number is always M = }/;mj;. Hence, the observer should not 
think that the gravitational matter are some extra particles, although the observer will notice 
its gravitational effect. It is obvious that our conclusions still hold when we consider the 
creation and annihilation of particles in the classical world. Because we use non-relativistic 
Schrédinger equation in this work to consider the classical approximation of gravitational 
effect, the application of Equation (31) should be confined to the scale of galaxies, which is 
sufficient for the present application. For larger scale of the universe, we should improve our 
model to consider the curved spacetime of general relativity. 

The fact that there are no unitary relations between Y’, at different times is not surprising. 
If we consider classical world branching, there are no unitary relations between Y; at different 
times even if there are no diffusion wave packets. Suppose we do not consider diffusion wave 
packets, then the only requirement for the existence of classical world at a given time t is 
< ¥.(t)|U(t, to) |¥u(to) >A 0, here ¥,,(to) is the universal wave function at an initial time 
to. For ¥,(t) this is similar, it is hard for us to answer why at time t we can construct such 
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Y(t) from ¥,(f) and ¥,,(t) and how to get the exact ¥’,(f), since we do not know the exact 
form of ¥/,,(f) at time f in the first place. However, so long as < ¥(t)|U(t, to)|¥u(to) >4 Ois 
satisfied, in principle we can regard ¥',(t) as a function of some significance. These analysis 
show that the sensible question to ask for future studies is what kind of initial universal wave 
function can cause ¥ ,(f) to be unorthogonal to ¥,,(f), where Y(t) contains ¥,(f). 

In the formal solution given by Everett [3], there are no longer any diffusion wave packets 
after some phase transitions occurred in the early universe in which particle wave packets are 
completely spread out. From the point of view of unitary evolutions, there is no such so called 
wave packet collapse, so the more natural hypothesis is that after phase transitions occurred 
there are mixtures of diffusion parts and correlated localized distribution of wave packets. The 
parameter d reflects the ratio of such mixtures. This inspires us that the process by which the 
universe formed its structures is a kind of phase transition in which diffusion and localized 
distributions coexist. After the phase transition takes place, during long periods of evolution 
of the universe the information about d is retained and it has two general characteristics: (i) 
since different classical worlds are formed from the same universal wave function, we expect 
that the values of d of Y., corresponding to other classical worlds should be about the same 
as our world; (ii) suppose at a time ft, we have d(t,) 4 0, then we think at a later time tz we 
have d(tz) = d(t,). Of course, this is not a proof, but an assumption of the universal value of 
the parameter d. It is similar to the formation of crystal that, the lattice constant of the crystal 
is always the same, although the initial condition to form the crystal can be significantly 
different, and we can not prove it by a calculation of the real evolutions. 


7. An Improved Schrédinger Equation for Classical World Wave Function 


The central purpose of the present work is to calculate evolutions of the classical in- 
formation {2j(t),--- ,Zn(t)} from unitary evolutions of Y,,(f), because only the classical 
information can be recorded by observers. The ideal way for this is to first solve the evolution 
of Y’,,(£), and search for the classical world at different times from Y(t). From the classical 
world wave function Y(t) obtained from Y(t), we can know the evolution of the classical 
information {X4(f),---,Ly(t)}. In the case of ¥p = Yy = Yu, we may get the classical 
approximation of the gravitational acceleration without using the method in this work. In 
the case of ¥. A Vy, Vy A Yu and ¥y(t1) F U (ty, t)¥,(t), however, special method should 
be used to consider the evolution of classical information {Xj(t),--- ,Zn(t)}. Among these 
three kinds of wave functions, only the universal wave function satisfies unitary evolutions of 
quantum mechanics. ¥,; and Y'y are only functions derived from Y,. The final result we get 
should focus on ¥, and its corresponding classical information. It is nonsense to calculate the 
long time unitary evolution of ¥y. 

In Figure 2 we summarize relationships between these three kinds of functions. There 
are unitary evolutions between the universal wave function at different times, we indicated 
this in the top of the figure for the set of universal wave functions {V1 (fo), Yu(t1), Yu(tz) } at 
times {fo,t1,t2}. From the universal wave function we can always search for classical worlds, 
here {Cx(to),Cx(t1),Cx(t2)} represents the set of classical information related by temporal 
causal relationships, while {‘¥* (to), ¥£(t1), YE (t2) } represents the set of classical world wave 
functions corresponding to {Cx(to),Cx(t1),Cx(t2)}. The so called classical theorem (here 
denoted by C) represents physical laws to predict how to evolve from Cx (to) to Cx (t1) and 
then to Cx (tz). Our analysis shows that, if we consider gravitational effects, then in general 
we cannot deduce good classical theorem from the unitary evolution of Y¢, since the best 
way is to derive {C,(fo),Cx(t1), Cx (t2)} directly from the set of universal wave functions 
{Vu(to), Pu(t1), Pu(t2)}. However, that is not a realistic approach so we introduced the set 
of functions {‘¥, (to), ¥(t), ¥x(t2)} to improve the classical theorem by considering more 
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carefully the gravitational interaction, compared with the classical theorem obtained from the 
unitary evolution of Y¢. 
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Figure 2. The relation between ¥,, Yc, ¥ and classical information. 


Observe that if we do not consider gravity, the solution of Y* from Schrodinger equation 
is a natural mathematical method to describe evolutions of the classical world, we only need 
to take proper account of classical world branching. Thus, we naturally consider the question 
of whether corrections from the x function would result in an improved Schrédinger equation. 
We define a new function Y* and let ¥*(t)) = YK(to). Then we think the following modified 
Schrédinger equation for ¥5(t) reflects evolutions of the classical world in a reasonably 
accurate manner: 


aye 





m A A _Pa(X,t 
= (% + Hesw + AGy = Gn yim; ) exes) i (35) 
j J 


Here 
R 1 x 
Hoy = —5 N 7; —?al’ 
1én=1 |"! 7 


M|Myn 


(36) 





while the last term on the right hand side of Equation (35) should be determined by astronom- 
ical observations. 

It is easy to show that the classical theorem derived from the above Schrédinger equation 
is the same as the classical theorem derived from combined consideration of ¥, and ¥,. We call 
Y5 gravity-rectified classical world wave function. Gy can also be regarded as a gravitational 
constant resulting from rectifying of gravitational constant of elementary particles. In the 
above equivalent Schrédinger equation, we do not have to restrict domains of variables 7; in 
Y5- The merit of the above improved Schrédinger equation lies in that almost all the physical 
phenomenon are fully considered, such as the classical world branching, all gravitational 
effects including the gravitational force due to diffusion wave packets, and quantum process 
in our laboratory. In a sense, we may regard this equation as the final result of our work. 
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So far we have introduced four kinds of functions: the universal wave function Y,(f), the 
classical world wave function Y(t), gravity-rectified classical world wave function Yo(f) and 
x function Y(t). Among them only the universal wave function is ultimate realistic entity 
and satisfies unitary evolutions of quantum mechanics. ¥,(f) is also realistic in some sense, 
since we can learn classical information of each particle in the classical world from ¥;,(t). 
In that sense Y,(t) is similar to ¥,(t), the difference being that when we solve Schrédinger 
equation restricted in the classical world domain UX;(t), Yo(t) gives better prediction of the 
evolution of the classical information than Y(t). From the point of view of perturbations, 
this is because Y,(t) takes better account of influences the universal wave function has on 
evolutions of the classical world. Without gravitational interaction, Y(t) and ¥;(t) are the 
same, and it is not necessary to consider Y(t) even when we still have ¥, A ¥. and ¥y A ¥y. 
As for ¥,(t), in some sense it is just a function without any realistic value, which is why we 
do not call it wave function; it is just a helper function to deal with gravitational effects using 
perturbation method and the diffusion part in Y(t) can change almost instantly. Hence, the 
classical world can not have memories of the information of {¢;(t), 6;(7;, t) } for each particle. 

Based on the above analysis, we see that evolutions of pq(7, f) in Equation (35) may not be 
calculated from unitary evolutions of Y(t). pg(7,t) may only be calculated from evolutions 
of the universal wave function. In practice, we need to observe the distribution of p,(7, t) 
through gravitational effects. Considering that Y(t) needs to be reconstructed at every 
instant, what we observe is actually the temporal and spatial average of o4(7,t), namely 











t+At/2 
FP a dt fd? ipa(y,t) 37) 
é VoAt : 
where Q in J q denotes the sphere with center 7, Vo denotes volume of the sphere. 
Thus the form of Equation (35) can be improved as: 
4 Oto (E) 3zPa(%,t) 
in = (Fo + Flan + By - Gu Eom ez ere 
(38) 


W(t) = U(t, to) ¥u(to) 
{¥5(t), Pult)} — pa(%, t) 


In the above equation the evolution of 64(X,t) over time is not determined by unitary 
evolutions of Y(t) (we have explained that Y(t) at different times cannot be related by 
unitary evolutions, not even short duration of time), but it is caused by the fact that the 
classical world with quantum correlation between nearby particles is contained within the 
universal wave function with diffusion wave packets all over the places. It seems plausible to 
assume that 6,4(x,t)/d is the mass density distribution calculated from Y(t). 


8. Conclusions 


As a summary, we discover that the universal wave function can be far richer than 
Everett’s many-world universal wave function. By analyzing the multi-body Schrodinger 
equation and structure of the universal wave function, we found that the universal wave 
function can have special and observable gravitational effects, and that makes the universal 
wave function a realistic entity. Everett’s formal solution is only the special model of our 
theory with d = 0. In this case, ¥; = Yy = Yu for a single classical world. As a comparison, 
most studies based on Everett’s idea are about the mechanism of classical world branching 
and decoherence based on the formal solution [3] , rather than the gravitational effect. 

Here we wish to emphasize that even though we provided a plausible solution to explain 
the gravitational effects of dark matter, our model does not exclude the possibility of existence 
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of new elementary particles [6] which also create extra gravitational contributions, such as 
axion model [11-13]. Based on quantum origin of the universe, we believe that the 7 function 
presented in this article are inevitable consequences of evolutions of the universe. Once 
those large scale diffusion wave packets and their gravitational effects in the universal wave 
function are verified, the existence of many worlds becomes a universal truth. Our researches 
provide new directions for future astronomical observation of the universe; that is, obtaining 
wave packets distribution information in the universal wave function through the observation 
of the gravitational effect of p,(7, f). 
Finally, we discuss several possible falsifiable predictions based on our theory. 


1. In the application of our theory to dark matter, we do not consider deviations of gravita- 
tional force due to fluctuations of yj and ¢;, and different 7; and ¢; for different particles. 
Because the gravitational constant is very small, all observed gravitational effects are 
for systems consisting of a huge number of particles. In this case, it is extremely hard to 
notice differences in gravitational effects caused by fluctuations of 7; and ¢;. It would 
be interesting to notice that after many years of the earth-based measurements [14] of 
Newtonian gravitational constant, there is an unexplained uncertainty in the measure- 
ment result of G [14]. One of the relevant experimental groups is asking for possible new 
physical principle [14] to explain the uncertainty. Whether this kind of fluctuations of 
y; and ¢; has chances to interpret the uncertainty of Newtonian gravitational constant 
is a question worthwhile of researching in the future. In future work, we will consider 
dedicated experimental scheme to detect this fluctuating gravitational effect, which has 
the chance to reveal the gravitational force due to diffusion wave packets, aside from 
astronomical observations. If the unexplained uncertainty in the earth-based measure- 
ments [14] of Newtonian gravitational constant can be found to have relevance to the 
fluctuation of 7y; and ¢;, it would be a strong support to our theory. 

2. In future work, we may consider the numerical calculations of Equation (38) to simulate 
the large scale structure of our universe, so that to give a reliable test of our theory. Of 
course, this needs the adjustment of the initial condition of the universal wave function, 
and has the chance to provides further clue for the initial condition of the cosmic inflation. 

3. | Our work not only shows that we need to change our understanding of gravitational 
effects in the classical world, but it also reveals that the Newtonian gravitational constant 
can be very different from the gravitational constant between elementary particles. 
Therefore, we also need to redefine magnitude of the gravitational constant between 
elementary particles, which provides the chance to test our theory in future. From the 
new gravitational constant G = (d + 1)Gy, we should reconsider the grand unified 
theory and may provide indirect method to test our theory. 

4. | We may consider an indirect test of our theory by simulating the model with ultracold 
atomic gases [15]. First we prepare a completely isolated quantum gas in which all parti- 
cle wave packets are diffusion wave packets. At the beginning we can even assume that 
there are no interactions between particles. Then, we artificially adjust the interaction 
energy to make phase transitions in this quantum gas. We use appropriate interaction 
energy to control the ratio of diffusion and localized parts after phase transitions took 
place. Under this condition, unitary evolutions can make the coexistence of different 
¥,(t) and ¥.(t) with definite structures. With rapid developments of techniques for 
ultra-cold atoms gas such as Feshbach resonance [15,16] to control the interaction be- 
tween atoms, we believe in the future there are prospects to simulate evolutions of the 
universe revealed in this article. 

When the temperature is extremely low, we may prepare the wave packets of ultra cold 
atoms gas with long range interactions contained within a potential [15] diffuse over a 
macroscopic distance. We can superpose a periodic potential on this system and position 
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all the atoms at bottom of the periodic potential [15,17]. Now we can employ interactions 
such as Feshbach resonance [15,16] to form some correlated definite structures (such 
as artificial crystal). Under this condition, the crystal structure is correlated with the 
diffusion parts. The interaction energy for this system becomes E = E, + Eq + E.g, where 
E.q is the interaction energy between correlated crystal and diffusion part, E, is internal 
interaction energy for the crystal and Ey stands for internal interaction energy of diffusion 
parts. Experimentally, people can measure the total interaction energy to determine 
if E.q really exists. Therefore this experimental scheme provides a way to simulate 
classical world wave functions presented in this article and the formation of classical 
world structures. Finally, even if we abandon motivations to simulate the universe, 
this problem per se is worthwhile of researching in the field of quantum multi-body 
manipulations. 


Note added: After the submission of our manuscript, the referee noticed a relevant paper 
[18] which argues that the solution of the Schrodinger equation including the gravitational 
potential may provide the possibility to explain dark matter. However, the mechanism to 
explain dark matter is completely different from our manuscript. In [18], the single-particle 
Schrodinger equation is used to explain the dark matter gravitational effect; while in our work 
the diffusion wave packet for the many-body universe wave function is shown to explain the 
gravitational effect of dark matter. 
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